Spin susceptibility of two-dimensional electrons in narrow AlAs quantum wells 
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We report measurements of the spin susceptibility in dilute two-dimensional electrons confined 
to a 45A wide AlAs quantum well. The electrons in this well occupy an out-of-plane conduction- 
band valley, rendering a system similar to two-dimensional electrons in Si-MOSFETs but with only 
one valley occupied. We observe an enhancement of the spin susceptibility over the band value 
that increases as the density is decreased, following closely the prediction of quantum Monte Carlo 
calculations and continuing at finite values through the metal-insulator transition. 

PACS numbers: 71.70.Ej, 73.43. Qt, 73.50.-h 
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As the density of a low-disorder two-dimensional elec- 
tron system (2DES) is lowered, the Coulomb energy dom- 
inates over the kinetic (Fermi) energy and is expected to 
lead to a host of correlated states at low temperatures. 
According to quantum Monte Carlo (QMC) calculations, 
for example, the spin susceptibility x oc g*m* of the 2DES 
should rise and diverge below a critical density at which 
the system attains a ferromagnetic ground state (5* and 
m* are the renormalized Lande g-factor and mass respec- 
tively) llj . Recent experiments [lllllllllllliliin 
a number of different 2DESs have indeed revealed a rise 
in g*m* with decreasing density, and there is even a re- 
port of the divergence of g*m* at a density where the 
system exhibits a metal-insulator transition Almost 
none of the measurements, however, agree quantitatively 
with the results of QMC calculations. Moreover, in Si- 
MOSFETs it is reported that the main contribution to 
the g*m* enhancement comes from an increase in m* 
rather than 17* also contrary to the QMC predic- 
tion. 

We report here measurements of g*m* and m* as a 
function of density in high mobility 2DESs confined to 
very narrow AlAs quantum wells (QWs). In several re- 
spects, this is a nearly ideal system to study the prop- 
erties of a dilute 2DES: it is a single- valley system with 
an isotropic in-plane Fermi contour and a thin electron 
layer. We find that g*m* in this system follows closely 
the prediction of recent QMC calculations . The agree- 
ment continues through the metal-insulator transition, 
and we observe no corresponding ferromagnetic instabil- 
ity at that density. Our measurements of m* exhibit a 
sample and cooldown dependence that betray a possible 
inadequacy of the Dingle analysis in this system. 

Bulk AlAs has an indirect band-gap with the 
conduction-band minima at the six equivalent X-points 
of the Brillouin zone. The Fermi surface consists of 
six, anisotropic, half-ellipsoids (three full-ellipsoids or val- 
leys) , with transverse and longitudinal band masses of mt 
~ 0.2 and m; ~ 1 respectively. This is similar to Si ex- 
cept that Si has six conduction-band valleys centered on 
points along the A-lines of the Brillouin zone. Electrons 
confined along the growth direction in a narrow (< 55 



A) AlAs QW occupy the single out-of-plane valley since 
its higher m* along z gives a lower confinement energy 
than for the in-plane valleys 0, . This is again 

similar to the 2DES in Si-MOSFETs fabricated on (001) 
Si substrates, except that in the Si case two valleys are 
occupied. 

We performed measurements on four samples (Al, 
A2, A3, Bl) from two different wafers (A and B). AU 
were Si-modulation doped, 45 A-wide, AlAs QWs with 
Alo.4Gao.6As barriers, grown on GaAs (001) substrates. 
We patterned the samples in a Hall bar configuration, 
and made ohmic contacts by depositing AuGeNi and al- 
loying in a reducing environment. Metallic front and back 
gates were deposited to control the carrier density, ra, 
which was determined from the frequency of Shubnikov- 
de Haas (SdH) oscillations and from the Hall resistance. 
Values of n in our samples were in the range of 0.54 to 
7.0 X IQ-'^^ cm~^, with a peak mobility of ^ = 5.1 m^/Vs 
in sample Al; this is comparable to the highest mobili- 
ties reported for Si-MOSFETs. Using the AlAs dielectric 
constant of 10, and the transverse band effective mass rub 
= 0.21 reported from cyclotron resonance measurements 
01, our density range corresponds to 2.8 < Ts < 9.6, 
where rs is the average inter-electron spacing measured 
in units of the effective Bohr radius. We made measure- 
ments down to T = 40 mK using low-frequency lock-in 
techniques. The samples were mounted on a tilting stage 
to allow the angle, 6, between the normal to the sample 
and the magnetic field to be varied in situ. 

In Fig. 1, we exhibit a characteristic trace of longitudi- 
nal resistivity at 6* = 0°. The high quality of the sample, 
as compared to previous 2DESs in narrow AlAs QWs 
[l^. is evident from the persistence of SdH oscillations 
to fields as low as 0.5 T. The oscillations at the lowest 
fields correspond to odd-integer Landau level (LL) filling 
factors (h') , as determined by the ratio of the Zeeman and 
cyclotron splittings which, in turn, is fixed by the prod- 
uct g* m* . A simple analysis reveals that the pxx minima 
at odd-integer i' will be strong at the lowest fields when 
4j+l < g*m* < 4j+3, for j a nonnegative integer; this is 
the case for our 2DES as we will discuss below. 

For a quantitative determination of g*m*, we measured 
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FIG. 1: Longitudinal resistivity {pxx) for sample Al at n 
= 5.55x10^^ cm~^. Inset: Temperature dependence of the 
amplitude {/^.p^x) of SdH oscillations, normalized to the zero- 
field resistivity (po), at v — 7 {B — 1.0 T) for samples Al (up 
triangles, n — 1.82x10^^ cm" 
triangles, n — 1.89x10^^ cm 

fits give effective masses of 0.45 and 0.30 respectively. 



''\ = 4.7 mVVs) and A3 (down 
/i = 3.3 mV Vs). The Dingle 



magnetoresistance in fields applied parallel to the 2DES 
plane {9 = 90°) to directly observe the total spin po- 
larization field at low densities (Fig. 2), and in tilted 
fields (0°<6'<90°) to determine ff*m* at low fields and 
high densities (Fig. 3). An example of our parallel field 
measurements is displayed in Fig. 2. The traces are sim- 
ilar to results from previous experiments 0, 0, 0, |Fj • a- 
metallic T-dependence is observed at zero magnetic field, 
followed by a field-induced transition to an insulating be- 
havior at high fields. There is also a resistance kink whose 
high field onset, marked by an arrow (Bp) in Fig. 2, sig- 
nals the total spin polarization of the carriers 0, 0, ^3 ■ 
In the inset of Fig. 2, we have displayed Bp versus n from 
our parallel-field measurements. By monitoring Bp as n 
is changed, we can extract g*m* as a function of n from 
the relation Bp — {2h/ e}n/ g*m* (ignoring the possibil- 
ity of nonlinear spin polarization, which we will address 
shortly); these g*m* are plotted in Fig. 4(a) (open sym- 
bols). The results are normalized by the band values, 
56=2 and m6=0.21 [ulEi- 

To augment our determination of g*m* at higher den- 
sities, we next performed measurements in tilted fields. 
The quantization of the 2DES energy into LLs depends 
on the component of the magnetic field perpendicular 
to the plane of the 2DES, while the Zeeman energy de- 
pends on the total field. By tilting the sample relative 
to the applied field, we can change the ratio of Zeeman 
and cyclotron splittings and therefore the strengths of 
alternating SdH minima The minima are weakest 

when a set of spin-down LLs cross a set of spin- up LLs. 
This corresponds to the condition, l/cos(0c) = 2i/g*'m*, 
where 9c is the angle at which the levels cross and i is the 



FIG. 2: Parallel-field magnetoresistance at a fixed density 
for different temperatures. A field-induced transition from 
metallic to insulating behavior occurs at 7.7 T, and the 2DES 
becomes fully spin polarized at Bp — 10.1 T (vertical arrow). 
Inset: Measured spin polarization fields (Bp) vs. density for 
samples Al (up triangles), A3 (down triangles), and Bl (cir- 
cles). The solid curve represents the quantum Monte Carlo 
(QMC) prediction for Bp and the dotted line the zero-field 
metal-insulator transition (MIT) observed in our samples. 



difference in LL index between the crossing levels. In Fig. 
3 we have plotted, at a fixed density, the strength of SdH 
minima, Ap^xi for several v as a function of l/cos(0); 
Apxx is defined as the average difference between the re- 
sistance at a SdH minimum and at the maxima on either 
side of it. If the maxima become indiscernible, such as 
when a minimum itself becomes a maximum, we simply 
take Apxx = 0. The LL crossings are indicated by the 
vertical dashed lines in Fig. 3, with the even-i/ minima 
experiencing a crossing at 9c = 57.7° and the odd-iy min- 
ima at 9c = 74.7°. The lowest i' at which LL crossings 
are observed is i^—2, giving i = 1 for the first observed 
crossing angle. The plot of l/cos(6'c) versus i, shown in 
Fig. 3 inset, confirms the above LL crossing condition. 
Using such analysis, we have extracted g*m*/ gtrnt versus 
n and plotted them in Fig. 4(a) (closed symbols). We 
note that the values of 17*771* derived from parallel and 
tilted field measurements are consistent with each other 
and with the condition discussed above for the odd filling 
SdH minima being strong at low fields. 

We now address the linearity of the spin polarization 
with field. In their study of dilute GaAs 2DESs, Zhu et 
al. and Tutuc et al. Q report a significantly nonlinear 
spin polarization as a function of magnetic field. Mea- 
surements by Tutuc et al. demonstrate that the finite 
thickness of the electron layer causes an enhancement 
of the single-particle m* in parallel field. This enhance- 
ment is significant when the magnetic length becomes 
comparable to or less than the 2DES thickness and ac- 
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FIG. 3: Strength of SdH minima as a function of tilt angle 
for several different filling factors. LL crossings are indicated 
by dashed lines, with the even-fillings experiencing a crossing 
at 8c = 57.7° and the odd at 74.7°. Inset: Dependence of the 
LL crossing angle, 9c, on crossing index, i. 



counts for the observed spin polarization nonlinearity in 
GaAs 2DESs in Ref. 0. In contrast, we do not expect 
significant finite-layer thickness effects for a 2DES in a 
narrow AlAs QW given its very small width. Our results 
indeed reveal no measurable nonlinear spin polarization 
with magnetic field. First, the LL crossings occur at the 
same 6c for all v of a given parity (Fig. 3) with only 
slight deviations at the very lowest i^. Second, the plots 
of l/cos(0c) versus i for various densities are all fit well 
by straight lines that pass through the origin (Fig. 3 in- 
set); this would not be the case if the rate of polarization 
changed significantly in field. Third, there is good agree- 
ment in our values of g*m* derived from tilted- field LL 
crossings and parallel-field Bp measurements, although 
the former measurement is performed at small, partial 
spin polarizations and low parallel fields while the latter 
is done at full polarization and high parallel fields. 

Having established the linearity of spin polarization in 
this system, we compare our Bp and g*m* with the pre- 
diction of a QMC calculation for an interacting 2DES 
which has zero thickness and is disorder free [ij. That 
prediction is shown in Fig. 2 inset and in Fig. 4(a) 
|17| . The agreement of the measured and predicted val- 
ues is surprisingly good, considering that the calculation 
ignores disorder. While we cannot rule out that such 
good agreement is fortuitous, we emphasize that there 
are essentially no adjustable parameters in this compar- 
ison. To convert the 2D density into rs, we used a band 
mass rrib = rut = 0.21; although there is some uncertainty 
regarding the precise value of the band mass in AlAs, the 
0.21 value is from cyclotron resonance measurements 0| 
and likely represents the most accurate value currently 
available. We add that we get good agreement between 
the QMC results and the measured 3*m* for 0.20 < mt 
< 0.23. This range is also consistent with a recent deter- 
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FIG. 4: Shown as a function of density are: (a) spin sus- 
ceptibility, (b) effective mass, and (c) mobility. In all plots, 
up triangles correspond to sample Al, squares to A2, down 
triangles to A3, and circles to Bl. The QMC prediction of 
g*m* at zero (dotted) and full (solid) polarization is included 
in (a), and the vertical dashed line indicates the density cor- 
responding to the zero-field MIT. 



mination of mt, from ballistic transport measurements in 
AlAs 2DESs [13 • On the other hand, SdH measurements 
of the m*, including our own, tend to give slightly larger 
values. As we discuss below, these latter measurements 
are problematic. 

In our samples, we observe a zero-field metal-insulator 
transition (MIT) at n=0.72 x 10" cm^^ (^^^8.3), indi- 
cated by the vertical lines in the inset of Fig. 2 and in 
Fig. 4. This is close to the value reported for wide AlAs 
QWs with comparable mobilities but with the in-plane 
valley occupied [191, and Si-MOSFETs with two out-of- 
plane valleys occupied ^3- More importantly, there is 
no apparent divergence of g*m* at the MIT in our sam- 
ples, and instead the measured g*m* coincide with the 
QMC curve continuously through the transition. The 
QMC calculation places the ferromagnetic transition in 
our system at n 1 x 10^" cm~^, well below the range 
accessible in our measurements. 

We also performed independent measurements of m* 
as a function of n [Fig. 4(b)] by monitoring the T- 
dependence of the amplitude of low field SdH oscillations 
and fitting the results with the Dingle formula (Fig. 1 
inset). Both the mass and the Dingle temperature were 
allowed to vary as fitting parameters. The masses that 
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we include in Fig. 4(b) at each n were obtained from 
data taken in the temperature range 0.3 to 0.8 K and 
fields in the range of 0.50 to 1.5 T. We have included 
two-sided error bars in Fig. 4(b) to reflect the spread 
of m* in this field range. Our measured m* at high n 
are comparable to those reported previously in narrow- 
well AlAs 0, At low densities, there appears to 
be some enhancement of m* that is similar in magnitude 
to that reported in Si-MOSFETs |1E3. however there 
is a significant spread in the m* values that depends on 
both sample and cooldown. As examples, in the inset of 
Fig. 1, we show Dingle plots for two different samples 
at nearly the same field and density that give markedly 
different values for m*. It is noteworthy that the discrep- 
ancies become large at approximately the same density 
where the 2DES mobility begins to drop rapidly [Fig. 
4(c)], though there is no direct correlation between the 
mobilities and m* values. The Dingle analysis may be 
producing anomalous results in this system, particularly 
as the density is reduced, for several reasons. At low 
densities, the zero-field resistance becomes increasingly 
T-dependent, complicating the application of the Din- 
gle analysis which presupposes no such dependence. We 
have tried the analysis employed in Ref. [Sj whereby the 
Dingle temperature is assumed to vary with T at the 
same rate as the zero-field resistance; this fails to appre- 
ciably reconcile the discrepancies between the m* values 
and simply reduces the overall enhancement at low n. 
The discrepancies could also be due to inhomogeneities 
in the 2DES density, which can vary between samples and 
cooldowns and which would become increasingly signifi- 
cant as the overall density is reduced. In any case, given 
that the discrepancies observed in m* are absent in g*m*, 
it appears that the values deduced from the Dingle anal- 
ysis of the SdH oscillations do not reflect the true m* at 
low densities in our samples. 

We emphasize that the values of 17*771* derived from 
parallel and tilted field measurements are consistent be- 
tween samples and cooldowns. This is despite differences 
in mobility, implying that the spin susceptibility is not 
very sensitive to the effects of disorder. The continued 
agreement of our measurements with the QMC calcu- 
lation at low densities, where screening of the disorder 
potential wanes and the role of disorder become more 
important [Fig. 4(c)], is consistent with this observation. 

We have argued that narrow AlAs QWs are ideal for 
measuring the role of interaction on the spin susceptibil- 
ity because of their single- valley occupation, thin electron 
layer, and isotropic in-plane Fermi contour. These same 
properties may be responsible for the differences between 
our results and those in other 2DESs. Measurements of 
g*m* in Si-MOSFETs generally find less enhancement 
over the band value than narrow AlAs QWs for most of 
the available density range. Initially, this is somewhat 
surprising considering that (001) Si-MOSFETs have a 
double-valley occupancy which should produce a more 



dilute 2DES. However, recent measurements in strain- 
tunable wide AlAs QWs find the same valley occupa- 
tion dependence of g*m* within a single material system, 
namely less enhancement in the two- valley case \2l\ . 

Another system of interest is the dilute GaAs 2DES 
where the electrons occupy a single, isotropic valley. In 
Ref. 1^, where g*m* was determined via Bp measure- 
ments, it was shown that the differences between the 
measured 17*771* and the QMC results could be essen- 
tially reconciled by taking into account the significant 
Fermi surface distortion that occurs in a thick electron 
layer system in the presence of a large parallel field. Zhu 
et al. 7], however, have reported that even in the limit of 
zero parallel field, the measured g*m* remains below the 
QMC results. While we do not have an explanation for 
this observation, we speculate that the weaker enhance- 
ment may stem from an effective reduction of the inter- 
action strength at zero field caused by the very large elec- 
tron layer thickness in the extremely dilute GaAs 2DESs 
of Ref. (2l] J rendering a system with a smaller effective 
than the one calculated assuming zero thickness. A full, 
quantitative understanding of the differences between the 
various material systems is still lacking. 
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